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We investigate the generation of gravitational waves in power-law inflationary models. The en- 
ergy spectrum of the gravitational waves is calculated using the method of continuous Bogoliubov 
coefHcients. We show that, by looking at the interval of frequencies between 10~^ and fO'' Hz and 
also at the GHz range, important information can be obtained, both about the inflationary period 
itself and about the thermalization regime between the end of inflation and the beginning of the 
radiation-dominated era. We thus deem the development of gravitational wave detectors, covering 
the MHz/GHz range of frequencies, to be an important task for the future. 
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I. INTRODUCTION 

Gravitational waves have the potential of providing us 
with a unique telescope and a unique source of informa- 
tion about the very early universe. 

Since the first calculation of the full spectrum of 
stochastic gravitational waves has been performed by 
Allen continuing and cornpleting work by other 
authors, including Grishchuk and Starobinsky 
Rubakov, Sazhin and Veryaskin [il], as well as by Fab- 
bri and Pollock 5], Abbott and Wise [1], and Abbott 
and Harari [?[, a lot of research has been done in this 
direction, particulary into exploring the potentialities of 
different inflationary models. Although the values ob- 
tained for the relative energy spectrum of the cosmologi- 
cal gravitational waves seem, at present, to be beyond the 
capabilities of the first generations of gravitational-wave 
detectors, this is a situation that is bound to be changed, 
given the effort that is being put into the planning of 
future detectors with considerably improved sensitivities 
like, for instances, NASA's Big Bang Observatory (BBO) 
Q and the Japanese Deci-hertz Interferometer Gravita- 
tional Wave Observatory (DECIGO) 0. 

An interesting question to ask is what can we learn 
about the physics of the very early universe, from the 
reading of the spectrum, assuming an ideal detector. In 
this article, we address this question with the help of a 
popular class of models, namely, power-law inflation. We 
shall show that, by looking at the interval of frequencies 
between 10~^ and 10^ Hz, and also at the GHz range, 
important information can be obtained, about both the 
inflationary period and the thermalization regime. Con- 
cerning the first range of frequencies, between 10^^ and 
10^ Hz, important work, along the lines of our paper, 
had already been done albeit by a different method, 
with similar conclusions to our own. Most of the de- 
velopments that have taken place in gravitational wave 
detectors cover the first range of frequencies, but work is 



also going on to study the possibility of constructing de- 
tectors covering the MHz and GHz range of frequencies 
[m . We deem these last developments to be very impor- 
tant, given the amount of information that can thereby 
be obtained, as we show with the help of a simple model 
for the transition between the inflationary and radiation- 
dominated eras. 

In this paper, we consider those gravitational waves 
generated by the usual mechanism of quantum fluctua- 
tions of the vacuum. In the MHz/GHz region another 
independent mechanism is possible, when the end of in- 
flation is followed by a period of parametric resonance 
(12i] . in which case the gravitational waves are directly 
sourced by the large inhomogeneities in the matter dis- 
tribution occurring during this period. In this paper we 
shall not address this important case, which has been 
investigated elsewhere [TslIT^. using both mechanisms. 

Throughout, gravitational wave production will be cal- 
culated using the method of continuously evolving Bo- 
goliubov coefficients. This method, applied to the pro- 
duction of particles in an expanding universe, was first 
investigated by Parker [l^. In the case of gravitational 
waves, the differential equations determining the contin- 
uous Bogoliubov coefficients were derived afterwards in 
a different, geometrical, way [l^. The method of con- 
tinuous Bogoliubov coefficients has advantages over the 
frequently used sudden transition approximation. Asso- 
ciated with the sudden transition there is always an over- 
production of gravitons of large frequencies, requiring an 
explicit cut-off for frequencies above the rate of expansion 
of the universe [l| . This cut-off is obtained in a natural 
way by the use of continuous Bogoliubov coefficients 
This is also a very practical method to calculate the full 
spectrum, from the very low frequencies corresponding to 
the present cosmological horizon, till those large, GHz, 
frequencies associated with the transition between infla- 
tion and the radiation-dominated universe. 

This paper is organized as follows. In Sec. |TT] we in- 
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troduce the equations defining the power-law inflation- 
ary model. The parameters of the inflationary potential 
are constrained using recent measurements from the cos- 
mic microwave background and large-scale structure. In 
Sec. mil reheating is incorporated in our model of evolu- 
tion through an elementary decay mechanism of the infla- 
tionary scalar field into a relativistic radiation fluid. De- 
spite its simplicity, our reheating potential may contain 
relevant features of more realistic potentials. In Sec. HVl 
we present the differential equations to determine the 
continuous Bogoliubov coefficients and address the rel- 
evant issue of the initial conditions to be used in the 
numerical integration of these equations. In Sec. |V] we 
describe our numerical simulations and present several 
gravitational-wave spectra, obtained for different values 
of the parameters of our model. The paper ends with the 
Conclusions. 



II. POWER-LAW INFLATION 

Let us assume that the evolution of the universe during 
the inflationary period is dominated by an homogeneous 
scalar field cj) with a potential 



1/(0) = Voe-^t 



(1) 



where Vq and A are positive constants. For a spatially 
flat Friedniann-Robertson- Walker metric. 



the Einstein equations for and 0(77) are 
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where a prime denotes a derivative with respect to con- 
formal time 77 and mp = l/\/G = 1.22 x 10^^ GeV is the 
Planck mass. 

The above set of differential equations admits the exact 
solution 



a{ri) = 01(771 - -qf, 



(6) 



where p < — 1, oi and 771 are arbitrary constants, and 



X = 4:\lTr rrip^ . 

P 



(8) 



Inflationary cosmology predicts a nearly scale- 
invariant power spectrum of density perturbations p^ . 



Within the slow-roll approximation, the power spectrum 
is given by |20|] 



Psik) = Psik,) ( A 



(9) 



where the scalar power-spectrum amplitude Ps{kc) and 
the spectral parameters ns{kc) and as(fcc) are evaluated 
at some pivot wave number kc- In this paper, we take for 
these quantities the values Ps{kc) = (2.45 ±0.23) x 10^^, 
ns{kc) = 1.0 ± 0.1 and |as(fcc)| < 0.04 at the pivot 
wave number kc = 0.05 Mpc~^, corresponding to dis- 
tance scales of the cosmic microwave background (CMB) 
and large-scale structure (LSS) [20| . 

The spectral parameters Ug and as can be defined 
in terms of the usual slow-roll parameters e, 77 andf , 
namely, Us = l-&e+2f] [2l| and = 16e?7-24e^-2^ d^. 
For the exponential potential H]), the slow-roll parame- 
ters are e = (p-|-l)/p, f] = 2{p+l)/p and ^ = 4(73+1)^/73^, 
yielding — ~{p + 2)/p and as = 0. 

Taking into account the above relation between Us and 
p and the CMB/LSS constraint on Us, we conclude that 
p ^ —2/1.9 « —1.053. Despite the fact that p takes 
values within a very narrow window, —1.053 ^ p < —1, 
it has quite a dramatic effect on the energy spectrum of 
the gravitational waves, as will be shown in Sec. |Vl 

Within the slow-roll approximation, the power spec- 
trum ([9]) can be written as 



Psik) 



1287r 



37n6 {dV/d(j)f 



P 



3773,^ p+1 



(10) 



where (/>c = </'('7c) is the value of the scalar field at the 
moment when the scale k^ exits the Hubble horizon dur- 
ing the inflationary period. Within our model, 0c can be 
chosen freely. Let us express it in terms of 4n = 4'{rii), 
the value of the scalar field at the time reheating begins. 
Taking into account that the number of e-foldings of ex- 
pansion between 0c and (pi is given by 



Nr. 



8tt f V 



Trip J dV/d(j) ' 



we obtain 



0c = 0j - Nc\j^^^m-p. 



(11) 



(12) 



Inserting this expression into Eq. ()10p . and using Eq. ([9|). 
we obtain a constraint on the infiationary scale Vq, 
namely, 

^0 = '-^Psik.) 



xexp<j4^/7r^-^-2iVc^j>m*, (13) 
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where 47 < A^c ^ 62 [20| and, as already mentioned 
above, 0i is a free parameter, marking the beginning of 
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FIG. 1: Potential 1^(0) for p = - 
2, implying A = 1.547m;7\ Vo 
8.152m" 



-1.05, <f)i — — 0.3mp and 
= 1.461 X 10-"mp, Uo = 



the reheating process. In the numerical simulations, de- 
scribed in detail in Sec. |Vl we will use (/i^ = — 0.3mp, 
Psikc) = 2.45 X 10"^ and = 55, implying that 
1/4^1^^ Tni6 foj. _i.o53 p < -1. 



Fq"'" 1.5 X IQis 



III. REHEATING 

In order to incorporate reheating in our model of evo- 
lution, we assume that the potential of the scalar field is 
no longer given by Eq. ([T|), but rather by 



Uo ie''^ - 1) 



4> ^4>i 
> 0i 



(14) 



where n is a even number and and v are chosen such 
that the potential and its first derivative are continuous 
at (/) = 0i (see Fig. [1]). 

This is clearly a toy potential. However, it may con- 
tain relevant features of more realistic potentials, arising 
from extensions of the standard model; these are likely 
to involve potentials defined by power-series expansions 
of scalar fields, or by exponentials [2^. An important 
example of the last type, may be seen in the application 
of the Salam-Sez gin six-dimensional supergravity model 
[23 | to cosmology [2^ . Our potential is similar to the one 
derived within the Salam-Sezgin model, if in Eq. we 
put n=2. 

Initially, the scalar field cj) has a large negative value. 
A period of pure power-law inflation takes place until cj) 
rolls down to values of the order of 0^. After that, the 
scalar field rapidly approaches the minimum of the po- 
tential, where V ~ L/q beginning to oscillate around 
it. The energy of the scalar field is then transferred to 
a relativistic radiation fluid and reheating takes place. 
The decay rate of the scalar field into the radiation fluid 
is governed by a phenomenological parameter, F^, which 
we choose to be of the order of 10~^mp — 10~®mp. 

During the period of transition between power-law in- 
flation and the radiation-dominated era, the evolution of 



the universe is described by the set of differential equa- 
tions 
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(15) 
(16) 
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(18) 



where prad is the energy density of the relativistic radia- 
tion fluid. 

The transition period ends when the energy density of 
the scalar held, p^ = 0'V(2a2) -|- ;7o(e'"^ - 1)", becomes 
much smaller than the energy density of radiation, prad- 
After that, we neglect the scalar field and the evolution 
of the universe till the present time is described by the 
set of differential equations 
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AttG 
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(19) 



(20) 



where w = 0.78 and the density of radiation, mat- 
ter and dark energy at the present time 770 are, re- 
spectively, prad,o = 4.6 X 10"^^ kg/m^, pmat,0 = 2.6 X 
10-27 kg/m^ and pdc,o = 6.9 x lO^^^ kg/m^ yielding 
R{^a) = 71 km s^^Mpc^^ 



IV. GRAVITATIONAL- WAVE SPECTRUM 

The dimensionless relative logarithmic energy spec- 
trum of the gravitational waves at the present time ryo 
is defined by 



f^Gw('7o) = 



1 



dpGW 



Pcrit(??o) rflnw 
3Jg'?(,o) "'^"°^^'^"°^' 



(21) 



where pcrit is the critical density of the universe, pcw 
and Lo are the energy density and angular frequency of 
the gravitational waves, respectively, iJ is the Hubble 
parameter and /3 is a Bogoliubov coefficient, such that 
|/3p = {Nk{r])) gives the number of gravitons. 

In order to calculate the amount of gravitons produced 
during the evolution of the universe we use Parker's 
method of continuous Bogoliubov coefficients [l^ [l^ . 
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according to which the Bogoliubov coefficients a{ri) and 
(3{ri), defined as continuous functions of time, are deter- 
mined by the set of differential equations 



2k 



'2k 



ae 



-2ik(ri — 7'ii 



a 
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(22) 
(23) 



where k — 2TTa{ri) / X[ri) — a{r])u;{ri) is the co-moving 
wave number. This method has been used by us in pre- 
vious works to investigate the generation of gravitational 
waves in different cosmological models [3, [iTi, i26il . 
Upon the redefinition 



the system of Eqs 

X" 



and (j23p become simply 



fc2 _ _ X = 



Y = -X'. 
k 



(24) 
(25) 

(26) 
(27) 



These differential equations can be solved exactly for 



pure power-law inflation. Indeed, for a{7]) oc (r/i 
Eq. ([26|) can be written as 



which admits the solution 



(28) 
(29) 



where y = fc(?7i — rf), ^ = ±(1/2 — p), and and are 
Bessel functions of the first and the second kind. 

The functions X and Y should satisfy the initial con- 
ditions \X{7] = — oo)| = 1 and \Y{ri = — oo)| = 1, corre- 
sponding to |q!(77 — ~oo)| — 1 and \l3{ri = — oo)| = 0. 
This implies that the integration constants should be 
ci — yJiT /2e''^ and C2 = ±i-\/7r/2e'*, where 5 is an ar- 
bitrary real constant. The solution of Eq. ([26]), which 
satisfies appropriate initial conditions, is then 



where H^^'^^ — 



X{y)^^^H\)'^^^{y)e\ (30) 
if2-p ~ Ji/2-p ± iYi/2-p denotes the Hankel 



functions. 

From Eq. ([27]) is now straightforward to obtain 



Y{y) 



(31) 



which satisfies the initial condition \Y{iri — — oo)| = 1. 



At the end of the pure power-law inflationary period, 
where our numerical integration starts {"q = rjj), the func- 
tions X and Y are given by 



(32) 
(33) 



where yi = k{rii — rji) = —kpa{rii)/a'{rii) and we have 
chosen the arbitrary constant 5 to be zero. Equations 
([5^ and ([55]) will be used as initial conditions for X and 
Y in our numerical simulations. 



V. NUMERICAL SIMULATIONS 

The numerical simulations begin at the end of the pure 
power-law inflationary period, rji, and continue up to 
the present epoch, 770. We use a Runge-Kutte method 
to solve the system of differential equations P^ - P^ . 
Equations p5)) and ([SO]) are used as constraint equations 
to check the accuracy of the numerical solution. 

As initial conditions we choose a{rii) = 1 and 4>{rii) = 
<f>i = — 0.3mp. Taking into account that, for ^ 0^, 
(^'V(2a2) = y(0)(p-|- l)/(2p- 1), we obtain 



^ 2(p+l) 
2p~ 1 



(34) 



where A and Vq are given by Eqs. ^ and p^ . respec- 
tively. Since any pre-existing radiation fluid is diluted 
during inflation, we choose p-cadiVi) — 0- Finally, a'{rji) 
is determined from the constraint (1151) to be 



(35) 



= y3vr^P,(fc,)e-^^(^'+i)/p 



where ip is the Planck time and we have taken into ac- 
count that t/o(e'"^- - 1)" = Voe-^'>''. 

Having determined the time evolution of a" /a from 
Eqs. ([Tl])-([lll), we then solve Eqs. ^ and ^ (again 
with a Runge-Kutte method) for different values of 
^{tIo) = k/a{'r]o) and with initial conditions (j32|) and 
(|55)) . Finally, we compute with Eq. (|25p and obtain 
^^Gw(?7o) from Eq. (^1]) . 

The gravitational-wave energy spectrum can be di- 
vided in three regions. In the first region, luq = 
Lu{rio) ranges from LOmin to Wrad^mat, where Wmin = 
27rc/dHubbic(?7o) ~ '^T^H{rio) = 1.4x10"^'' rad/s is the an- 
gular frequency of a gravitational wave with a wavelength 
equal, today, to the Hubble distance, and Wrad^mat is to- 
day's value of the angular frequency of a gravitational 
wave which had a wavelength equal to the Hubble dis- 
tance at the time when the energy density of radiation 
became equal to the energy density of matter (corre- 
sponding to the transition from a radiation-dominated 
to a matter-dominated universe). Because of red-shifting. 
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FIG. 2: Time evolution of ff" during the transition from the 
inflationary period to the radiation-dominated era, for ojo = 
lO" rad/s, p = -1.05, = -0.3mp, = 5.0 x lO'^mp and 
n = 2. 



FIG. 3: Time evolution of during the matter- and dark 
energy-dominated eras, for ljo = 5 x 10~^^ rad/s, p = —1.05, 
(/>i = — O.Smp, = 5.0 X lO^'^mp and n = 2. 



^lad^mat IS of the Order of 10 rad/s. In the second (in- 
termediate) region, loq ranges from Wrad^mat to w^^rad, 
where w^^rad is today's value of the angular frequency 
of a gravitational wave which had a wavelength equal to 
the Hubble distance at the time when the energy den- 
sity of the inflaton became equal to the energy density 
of radiation (marking the end of the transition period 
between the inflationary and radiation-dominated eras); 
this frequency is of the order of 10^ rad/s. Finally, in 
the third region, ujq ranges from cj^^^iad to (ximax, where 
the maximum angular frequency, which is of the order of 
10"'^° rad/s, is today's value of the angular frequency of a 
gravitational wave which had a wavelength equal to the 
Hubble distance at the end of the inflationary period. 

Let us now show, using the method of continuous Bo- 
goliubov coefflcients, that in the second (intermediate) 
region of the spectrum, the relative logarithmic energy 
density of the gravitational waves, contrarily to the situ- 
ation in exponential inflation, depends on the frequency 
of the waves, namely, f^cw oc cjg^^'*'^'' . 

Taking into account that, according to our numerical 
simulations, the scale factor at the present time is of the 
order of lO'^" [for a(ryi) = 1], and using Eq. ([55)1 to eva- 
luate a'{rii), we obtain 



yi « -3.5 X lQ-^"p ' 



p+1 



.e55(p+i)/p ujo s. 



(36) 



Let us consider values of p in the interval —1.053 ^ p ^ 
-1 - e. From Eq. follows that yi < I for ujq ^ 2.2 x 
10^ rad/s and e = 1.9x10^"', allowing us to expand X{rji) 
and Y{r]i), given by Eqs. ([5^ and ((551) . i^i converging 
power series. Keeping just the leading term, X{rji) oc 

and Y{Tf) oc yi~^, we obtain that 0^{rii) <x. yf'^ cx 

, 2(P-1) 

In order to determine the subsequent time evolution 
of we have to solve Eq. (pS)) . Let us consider sep- 
arately three periods: i) the transition from inflation 
to the radiation-dominated era {jji rj < ri,f,^ii^d) ii) 



the radiation-dominated era (jy^^rad ^ V < '7rad^mat) 
and Hi) the matter- and dark energy-dominated eras 

(?7rad->mat ^ V ^ Vo) ■ 

During the transition from inflation to the radiation- 
dominated era, a" /a is a rather complicated function of 
time and, therefore, Eq. ((26)) must be solved numerically. 
In Fig. [5] the time evolution of is shown for ojo = 
10^ rad/s, p = -1.05, 4>i = -0.3top, F^ = 5.0 x lO^^mp 
and n = 2. For other values of luq, the time evolution of 
is the same, provided that ujq ^ 2.2 x 10*" rad/s. In- 
deed, for such values of ujq, P = a^oj^ < 1.4 x 10^^^ ip^ 
is always smaller than a" {rji) / a{rji) . To see this, we in- 
sert 4>'{rii) from Eq. ((34)) into Eq. ((T6)) . use the fact that 
C/o(e'"^' - 1)" = Vbe"^"*', with A and Vo given by Eqs. (0 
and ((T5)l . obtaining, for e = 1.9 x 10"*^, 



— {r]i) = 37r 
a p{2p-l) 



Ps{kc)e- 



> 2.9 X IQ-^^tp^. 



(37) 



In the first stages of the transition between the inflation- 
ary and radiation-dominated eras, a" /a increases and the 
relation fc^ < a" /a continues to hold. Therefore, we can 
neglect in Eq. ()26p . thus obtaining the same time evo- 
lution of for different values of ujq. To the end of the 
transition, a" /a approaches zero and graviton production 
ceases, with becoming constant (see Fig. [2]). 

During the radiation-dominated era, remains con- 
stant. Indeed, during this era the scale factor is propor- 
tional to the conformal time, implying a" /a = 0. Then 
Eqs. dMl) and admit the solution 



X{r]) = X(?70^rad) cos 9(77) - «y(?70^rad) singf(?7), (38) 
Y{ri) = y(770^i.ad) cos 9(77) - iX(770^i.ad) singf(?7), (39) 



where 5(77) = k{ri — ?70-irad)- Now, using Eq. ((^5)) . it is 
straightforward to show that /3^(?7) — /3^(770^rad)- 

During the matter- and dark energy-dominated eras, 
remains constant, if luq ^ 10~^^ rad/s. Indeed, for 
such values of ljq, our numerical simulations show that 
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FIG. 4: Gravitational-wave spectra for different values of p: 
p — —1.01 (violet, upper curve at high frequencies), p = —1.03 
(blue, middle curve) and p = —1.05 (red, lower curve). In all 
spectra we have used (f)i — — 0.3mp, — 5.0 x lO^^mp and 
n = 2. 



is always greater than a" /a and, once again, we can 
neglect a" /a in Eq. ([^5]) . obtaining, similarly to the 
radiation-dominated era, that remains constant, i.e., 

=/3'(^/rad^mat). 

Putting everything together, we arrive at the conclu- 
sion that, for 10^"'^^ rad/s ^ wq ^ 10^ rad/s (interme- 
diate region of the gravitational- wave spectrum), gravi- 
tons are produced during the transition from the in- 
flationary period to the radiation-dominated era, but 
not during the radiation-, matter-, and dark energy- 
dominated eras, and that this production does not de- 
pend on the value of wq- This, together with the fact that 
/3^(r/i) cx ujI^^~^\ implies that /3^(?7o) « '^o^''""^^- Now, 
using Eq. (j2ip . we finally obtain that flowivo) Wq'^~'^^\ 
i.e., for 10~^^ rad/s S S rad/s the relative loga- 
rithmic energy density of the gravitational waves depends 
on the frequency of the waves, in agreement with the re- 
sults obtained by Sahni using a different method p^ . 
This contrasts with the results obtained for exponential 
inflation, where the spectrum in the intermediate region 
is flat. 

Let us point out that for 1.4 x 10~^^ rad/s S ^t) S 
10~^^ rad/s, the energy spectrum is steeper than in the 
intermediate-frequency region. This is due to the fact 
that, in the beginning of the matter-dominated era, a" /a 
is greater than fc^ for such range of frequencies, leading 
to an extra production of gravitons. After a while, a" /a 
approaches zero and graviton production ceases, with f3^ 
becoming constant (see Fig. [3]). 

We now produce a few gravitational-wave spectra, 
coming from our numerical simulations. For each spec- 
trum, we have to specify the values of p, (pi, n and F^. 
While p is constrained by measurements of the cosmic 
microwave background and large-scale structure to lie in 
the interval —1.053 ^ p < —1, the latter three parame- 
ters can be chosen freely. Our choice of p fixes Vq and A 
and, thus, the level of the spectrum. 
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FIG. 5: Gravitational-wave spectra for different values of n: 
n = 2 (red, lower curve), n = 4 (blue, second curve from 
below), n = 6 (violet, third curve from below), n — S (orange, 
second curve from above) and n = 12 (black, upper curve). 
In all spectra we have used p — —1.03, = — 0.3mp and 
= 5.0 X 10"*mp. 
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FIG. 6: Zoom of the high-frequency region of the spectra 
shown in figure [5] 



Gravitational-wave energy spectra, for different values 
of p, are shown in Fig. [H For intermediate values of 
the gravitational-wave frequency (10""'^^ rad/s ^ ^ 
10'' rad/s) the spectrum has a slope that depends on the 
value of p. When p approaches —1, the model becomes 
closer to exponential inflation and the spectrum closer 
to flat, as expected. For luq < 10"-"^^ rad/s, the energy 
spectrum is steeper than in the intermediate-frequency 
region, as explained above. The behavior in the high- 
frequency region of the spectra {ojq > lO*" rad/s) can be 
understood as follows. Near its minimum, the potential 
V{(l)) is proportional to 0". An homogeneous scalar field, 
with a potential V{4>) oc 0", oscillating rapidly relatively 
to the expansion rate of the universe, behaves as a fluid 
with an equation of state p = jp, 7 depending on n in 
such a way that the fluid becomes more and more stiff 
as n increases [131 • Furthermore, the relative logarith- 
mic energy density of the gravitational waves decreases 
at high-frequencies if inflation is followed by a period 
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dominated by a fluid with a dust-like equation of state 
(p — 0), while an increase of the relative energy den- 
sity is observed at high frequencies if inflation is followed 
by a stiff matter phase {p = p) [2^ . Since the spectra in 
Fig. [4] were obtained for n = 2, in which case the rapidly- 
oscillating scalar field behaves like dust, a decrease at 
high frequencies of J^gw is expected. 

The dependence on the parameter n of the relative log- 
arithmic energy density of the gravitational waves at high 
frequencies is clearly illustrated in Figs. [5] and [51 where 
the gravitational- wave energy spectra for different values 
of n are shown. The part of the spectra corresponding 
to frequencies between 10'' rad/s and 10^" rad/s does re- 
flect some of the characteristics of the transition between 
the inflationary and radiation-dominated eras, namely, 
an increase in the value of the parameter n corresponds 
to an increase of the relative logarithmic energy density 
of the gravitational waves produced during the period of 
thermalization near the minimum of the potential. 

Measurements of the cosmic microwave background 
radiation impose, for angular frequencies corresponding 
to the present size of the cosmological horizon Whor = 
1.4 X 10~^^ rad/s, an upper limit on the relative loga- 
rithmic energy density of gravitational waves, namely, 
f^Gw(whor, r?o) < 1-4 X 10"^° 0. All the spectra pre- 
sented in this paper satisfy, by far, this condition. In 
fact, one could even relax the constraints on p and 
Vb, discussed in Sec. |TT1 and still obtain gravitational- 
wave energy spectra satisfying the above constraint on 

VI. CONCLUSIONS 

In this work, we have investigated the generation of 
gravitational waves in power-law inflationary models. 
The parameters of the inflationary potential were con- 
strained using recent measurements from the cosmic mi- 
crowave background and large-scale structure. We have 
incorporated reheating in our model through an elemen- 
tary decay mechanism of the inflationary scalar field into 



a relativistic radiation fluid. Despite its simplicity, the 
reheating potential may contain relevant features of 
more realistic potentials, as was emphasized through an 
approximate comparison with the potential derived in 
Ref. .^M from the Salam-Sezgin model. 

We have used the method of continuous Bogoliubov 
coefficients to calculated the gravitational-wave energy 
spectrum for different values of the parameters of our 
model. All the spectra we have obtained satisfy the 
constraint imposed by measurements of the cosmic mi- 
crowave background radiation at angular frequencies cor- 
responding to the present size of the cosmological hori- 
zon. Between 10~^^ rad/s and 10^ rad/s the spectrum 
is not ffat, but has a slope that depends on the value of 
p, and so on the value of the parameter A defining the 
inflationary potential (Fig. 2]). When p approaches —1, 
the model becomes closer to exponential inflation and 
the spectrum closer to flat, as expected. The influence of 
the inflationary regime does extend a long way towards 
low frequencies and the observation, or non-observation, 
of such a slope could, in principle, give us information 
on the type of inflation. On the other hand, the part of 
the spectrum corresponding to the frequencies between 
10^ rad/s and 10^" rad/s does reflect some of the char- 
acteristics of the transition between the inflationary and 
radiation-dominated eras, in particular the behavior dur- 
ing the period of thermalization near the minimum of the 
potential, as explained at the end of the preceding sec- 
tion, and shown in Figs. [5] and [6l This is why we believe 
it is important to continue work already done to develop 
detectors covering the MHz and GHz range of frequen- 
cies: they might give us information difficult to obtain by 
any other means. This was one point we wanted to make 
with our present work, notwithstanding the fact that the 
model we used was a very simple one. 
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